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Abstract 

We are presenting an algorithm capable of simplifying tensor polyno- 
mials with indices when the building tensors have index symmetry prop- 
erties. These properties include simple symmetry, cyclicity and those due 
to the presence of partial and covariant derivatives. 

We are also including some examples using the Riemann tensor as a 
paradigm. 

The algorithm is part of a Mathematica package called Tools of Tensor 
Calculus (TTC)[web address: http://baldufa.upc.es/ttc 



1 Introduction 

The two main languages commonly used in writing tensor calculus expressions 
are the intrinsic notation and the index notation. The intrinsic notation seems to 
be the preferred language for making computer algorithms work with symbolic 
tensor expressions. The index notation, on the other hand, is extremely powerful 
in expressing and manipulating tensors. There are some simple expressions (in 
index notation) which are difficult to express in intrinsic notation for the sake 
of introducing new operators with new properties. This is the case such as 
expressing the simple operation of raising an index like (|f|) 

Ti j k =g^T imk (1) 

(summation over repeated index is understood). 

Recently some computer tools capable of working with index notation at a 
symbolic level have been appeared [|] |[ |j. 

fn the present study we are facing to the problem of tensor polynomial sim- 
plification using index notation, fn a recent article H we dealt with this problem, 
restricted to the case of simple index symmetry properties of the building tensors 
like 

Ti„ a) ...i„ M + aT h ... in = (2) 
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a being a scalar and a a permutation of {1, n\. (||) includes symmetry and 
antisymmetry of any sequence of indices, and pairsymmetry. We will refer to 
these kind of properties as homogeneous monoterm. The algorithm proposed 
in pi does not use any kind of library, neither temporal nor permanent and 
for that reason it does not require a lot of computer memory. V.A. Ilyn and 
A.P.Kryukov || propose a method capable of working with index symmetry 
properties which includes cyclic type properties like 

Tijk + Tkij + Tjki = (3) 

These were included by the authors in a set called multiterm like properties. 
Although the problem was satisfactorly solved theoretically, the computer time 
and the amount of computer memory needed to make elementary simplifications 
was large enough to consider that his study may have some practical and fun- 
damental limitations ( since an expression with 1 1 index needs about 2500Mb) . 

One of the problems encountered by V.A. Ilyn and A.P.Kryukov is related to 
the dummy indices. They consider the possibility of renaming dummy indices 
as new relations to add to the set of relations generated using the true index 
properties of building tensors. So in the monomial as SjkT^ 1 as well as the 
properties due to the fact that S is a symmetric tensor they add the ones 
deduced from the fact that the indices j, k are dummy, that is to say 

SjkT^ - S kj P ki = (4) 

The number of this kind of relations increase factorially with the number of 
dummy indices, so the algorithm needs to work very hard for monomials of 
interest. 

The algorithm that we present here, written in Mathematical language, is 
a part of the package Tools of Tensor Calculus (TTC)M, ^ || It can work 
with monomials whose building tensors have both monoterm and multiterm 
properties. In addition to dealing with dummy indices it also handles symmetry 
(including the partial derivatives) and antisymmetry of any sequence of indices, 
cyclic type properties and the Ricci properties due to the presence of two or 
more covariant derivatives. 

Taking the Riemann tensor as a paradigm our algorithm can handle mono- 
mials that have this tensor as a factor with the known properties: 



Rijkl = —Rjikl = —Rijlk = Rklij 
Rijkl + Riljk + Riklj — 

Rijkl, m ,n — Rijkl, n ,m 

Rijkl;m Rijmk;l ~t~ Rijlm;k 

Rijkl ;m....;n = aRijkt-,; n....; m + (Riemann terms) (5) 
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We have added the possibility of introducing simplifying properties such as 

Rm imj = R ij i Rm m = R ( 6 ) 

explaining how to do that. The above properties are very usefull because, al- 
though they may not be needed, from the theoretical point of view, they allow 
drastic reduction the amount of indices used by the algorithm , during calcula- 
tion. 

We have left the inclusion of the properties derived from the dimension of 
the base space, for the Riemann tensor and for the generic case for a future 
study. 

2 The problem 

The problem can be stated as follows 

Given a tensor monomial written in index notation ( dummy index included) 
together with all index properties of the building tensors , find a unique canonical 
equivalent polynomial. 

The canonical equivalent polynomial is derived from the complete set of 
monomial relations generated from the exhaustive application of all index prop- 
erties of the tensor factors over the original monomial. The solving procedure 
is related to some ordering criterion of the set of monomials. 

Dummy indices 

Let Mq(If, Jn) be a given tensor monomial, with F free indices Ip = ...,ip} 
and D dummy indices Jg = {j a (i), j a (D)} , a being any permutation of 
{1, ...,£} 
Example: 

M ({i, k, I}, {j, m,p}) = S ijk T{™T p lm (7) 

Due to the meaning of the dummy indices, if I is the identity permutation 
the following identities hold 

M (If,J%) = M (If,J I d ) Va (8) 

This means that there will be many equivalent expressions for a given monomial. 
This is the first problem to be solved. 

In || we looked at a method in which the monomial could be rewriten in a 
canonical form with respect to the dummy indices independently of the num- 
ber of dummy indices present in the monomial, without using many computer 
memory. That is to say, from all the lists only one is necessary, Jp>. 
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Essentially the method consists of using an ordering criterion that can be 
applied to a monomial without prior consideration of all the monomials to be 
ordered. 

We will call PCanonic [0] the function implemented for the algorithm. We 
have 

M (I F ,J%) ?C ^ icm M (I F ,J D ) (9) 

In the following sections we will handle monomials in the form Mo (If, Jd)- In 
this way we will not mention the dummy indices problem anymore. 

Tensor properties 

Of all tensor properties we will emphasize those that will be of importance 
within the algorithm. We are referring to the properties of single tensors of a 
given monomial. 

• Monoterm properties: the result of their application over a given mono- 
mial is also a monomial. 

• Multiterm properties: the result of their application over a given mono- 
mial is a polynomial (more than one term). 

• Homogeneous or non-ordered properties: they are properties of only 
one tensor whose application over a given monomial does not guarantee a 
better (alphabetical) ordering result without considering other properties. 

Example 

If we have Sij — Sji = or Tijk + Tkij + Tjki = in principle we do not 
know how to use it to improve the ordering of the monomial SijkT^T™ . 
We must take all the relations generated from the above-mentioned prop- 
erties into account in order to know which of the resulting monomials or 
polynomial is the best ordered. 

• Inhomogeneous or ordered properties : These arc properties relating 
more than one tensor, which we know how to use to improve the order of 
a given monomial. 

Example 

The property T™ i — Ai guarantees an improvement of the ordering of the 
monomial S ljk TP pm T™ 

These kind of properties can be used directly over the working monomial. 

Not all the property falls into these last two definitions. Clear counterex- 
amples are the Ricci relation 

V t . m . n - Vi. n . m = (]/(Riemann terms)) (10) 
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From (|TC]) we see that the right hand side terms can always be considered 
more ordered than the left hand side terms, although we do not know 
which of the left terms is the most ordered. In this case we must take into 
account all the possible substitutions of the homogeneous part in order to 
generate the properties of a given monomial. 

The consideration of an ordered property can make properties that in prin- 
ciple are not ordered become ordered or partially ordered properties. A 
clear example is the definition of the Ricci tensor in terms of the Riemann 
tensor R m im j = Rij that turns the character of the Bianchi II identity of 
the Riemann tensor resulting 

K"ijk;m = Rki;j ~ Rij;k (11) 

a partially ordered property, since we can consider that the terms con- 
taining the Ricci tensor are more ordered than the rest. If we want to 
include the definition of the Ricci tensor in the set of properties of the 
Riemann tensor we must include ([Tl|), as well as the properties belonging 
to the Ricci tensor itself. Although such a procedure can be difficult and 
unnecessary because we can work exclusively with the Riemann tensor, it 
is very convenient to eliminate superfluous indices. 



How to handle monoterm properties 

Given a monomial Mq{I f , Jd) the exhaustive application of the monoterm prop- 
erties generate a set of relations of the type 

M (I F , Jd) + a oi Mi(I F , J D ) = 
Mi(I F , J D ) + aijMjilp, J D ) = i,j = 0, 1, ... (12) 

where a i, ay are scalars. This set of relations can always be solved in favour of 
the most ordered monomial. The result will be 

Mi(J F , Jd) = biM s (I F , J D ) V^ = 0, 1, ... (13) 

hi being scalars and M s (I F , Jd) the most ordered monomial. The specific way 
to solve ( |T2| ) will be explained in section(3) . 

We will call PCanonic [1] the function implemented by the algorithm 

Mi(I F , J D ) PCan =5S> ic [1] biM s {I F , J D ) V* = 0, 1, ... (14) 
In the following section we will handle monomials in the form M S (I F , Jd). 
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How to handle multiterm properties 

Given a monomial Mq (If, Jd) the exhaustive application of the multiterm 
properties generates a set of relations of the type 

M*{I F , Jd)+Yj t>oiM?(I F , J D ) = 

i 

M?{I F ,J D ) + Y,h j M?(I F ,J D ) = i,j = 0,1,... (15) 

3 

where b i, bij are scalars. 

This set of relations can always be solved in favour of the set of the most 
ordered monomials that we will denote by Mf c (I F , J D ). In this way given any 
of the monomials Mf (I F , Jd), generated from Mq (I F , Jd) by the multiterm 
properties, we have 

Mf {Ip, J d ) = ^ C H M j° i 1 ? ,Jd) (16) 

3 

Cij being scalars. 

The specific way to solve (15) will be explained in section(3) . 

We will call PCanonic [2] the function implemented by the algorithm 

M?(I F ,J D ) PCan 3 icC2] ^ C4 ,M/ c (/ F ,J I5 )V l = 0,l,... (17) 



3 Algorithm description 

The algorithm for simplification of polynomial tensors in index notation, called 
Simplif yAlllndex, has the following functions: 



Functions acting on single tensors 

• InputTensor: 

InputTensor [Tsymbol, basis, ranklist] declares that the symbol TSymbol 
will be a tensor in the basis basis and with rank and type ranklist. 

Example 

InputTensor [T , XX , {1 , 1 , 1 , 1>] 

• InputSymmetries [Tsymbollindex] ,symmetrispecijications~\ declares that 
the tensor, previously introduced using InputTensor, Tsymbol has the 
symmetry properties of the indices, positioned as indicated by the argu- 
ment index, according to the specifications symmetriespecifications. 

Example 1: 
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InputSymmetries [T [i , j ,k, 1] , {i ,k, 1}] declares that T is symmetric 
with respect the index {i,k,l} 

Example 2: 

here we declare the same properties as the Riemann tensor, see (^|), for 
the tensor T 

InputSymmetries [T [i , j ,k, 1] , 

{{i, j}},{{k,l}},-C-Ci,j>,{k,l}}, Cyclic [j , k, 1] ] ; 
InputSymmetries [T [i , j ,k, 1, . ;m] , Cyclic [k, 1 ,m] ] 

• BasicRules [n] : once the symmetries of indices have been declared the 
corresponding rules that implement these symmetries are stored in the list 
BasicRules [n] . If the property is monoterm n—1 and if it is multiterm 
n—2. The value n=0 is reserved for the case of ordered properties which 
are generally introduced explicitly. 

Functions acting on monomials 

In the following we will assume that there are no ordered properties stored in 
BasicRules [0] as they are unnessencial. 
The main function is Simplif yAlllndex: 

Index [me£ricname,Simplif yAlllndex] [polytens] simplify, i.e. canonize, 
the polynomial polytens taking into account all the declared properties of the 
building tensor, as well as the dummy indices and the properties related to par- 
tial and covariant derivatives. Simplif yAlllndex is applied over each term of 
polytens. This function has, in turn, the following internal functions (We will 
suppose that the initial monomial is M (Ip, J£)): 

• PCanonic [0] : first we canonize the monomial with respect to the dummy 
indices. This step is taken each time a monomial is the result of some kind 
of index manipulation. So our initial monomial and all other generated 
monomials will have the form Mi(Ip, Jo) 

• LlistaMonomis [n~\ [count [n] ] : in these lists all the monomials considered 
before the present case are stored . n=l for monoterm properties and 
n=2 for multiterm properties. count[n\ is a counter which says that the 
monomials therein are the count[n~\ -th studied monomials. 

• LlistaRules [n] \_count\_n~\~\: in these lists the simplifying rules of the 
monomial in LlistaMonomis [n] [count In] ] are stored. 

• PCanonic [1] canonize a monomial with respect to the monoterm proper- 
ties. First look at the lists LlistaMonomis [1] [cownt[l]]. If the working 
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monomial is in one of these lists it applies the corresponding rule stored 
in LlistaRules [1] Lcountllll . If this is not the case it applies the rules 
stored in BasicRules [1] exhaustively generating a list of new properties 
which will be solved by PFind giving as a result LlistaMonomis [1] [countLl] ] 
and LlistaRules [1] [count[l~\ ] and it applies these rules to the working 
monomial. 

• PCanonic[2] canonizes a monomial with respect to the multiterm prop- 
erties once it has been canonized by PCanonic [1] . First look at the lists 
LlistaMonomis [2] \_count[2~\ ] . If the working monomial is in one of these 
lists it applies the corresponding rule stored in LlistaRules [2] [count[2~\ ] . 
If this is not the case it applies the rules stored in BasicRules [2] exhaus- 
tively generating a list of new properties which will be solved by PFind giv- 
ing as a result LlistaMonomis [2] Lcountl2l1 and LlistaRules [2] [cotmt[2]] 
and applying these rules to the working monomial. 

• PFind is a simple tool that solves any system of linear equations with a 
fixed ordering criterion. 

Let {x\, X2 1 ■ xm} be the space of ordered variables to be solved. Let 

eqlist = {c 1 +a y x J = 0,c 2 + a 2j Xj = 0, c N + a Nj Xj = 0} j = 1,...M 

(18) 

be the list of N equations with M variables ordered first taking into ac- 
count the variables and then the coefficients. That is to say, if in the n-th 
equation the most disordered variable is x rn then x m+ \ cannot be present 
in the (n — l)-th equation. 

PFind takes the initial list of equations and starts the following loop: 

Step 1: PFind takes the first equation in eqlist. It is solved with re- 
spect to the most disordered variable. This result is stored in eqrule [1] . 
eqrule[l] is applied over eqlist and after eliminating zeros this result 
is assigned to eqlist. Go to step 2 

Step n > 1: If eqlist is empty the loop is finished. On the contrary 
PFind takes the first equation in eqlist. It is solved with respect to the 
most disordered variable. This result is stored in eqrule [ri] . eqrule [n] 
is applied over eqrule [i] i = l,...,n— 1. eqrule In] is applied over 
eqlist, the zeros are eliminated and the result assigned to eqlist. Go 
to the step n + 1. 

eqlist will remain empty if the system is compatible. These must be the 
case for a set of compatible tensor properties. 

The computer time and memory of this algorithm have been controlled. The 
session presented in the appendix can be performed with no more than 2.6Mb. 
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4 Conclusion 



We have solved the problem of simplification of tensor polynomials in index 
notation with respect to a large set of properties for the building tensor satis- 
factorily. The computer time and memory needed are sufficiently reasonable so 
we think that the program can be of interest for practical purposes. 

It is not our intention to say that TTC is a great program. Within Math- 
ematica there are very good programs (See MathTensor |8j ) Nevertheless we 
think that in the field of simplifying tensor polynomials we have reached a high 
level performance. Specifically, due the generic character of the algorithm pre- 
sented (the only requirements is that the input polynomial follow the usual 
mathematical rules) makes that it can be useful in the kind of fields where ten- 
sors are a natural language. As is well know tensor calculus is power enough to 
cover large fields of mathematics, physics and engineering. 

There is an important set of properties not included in the program. These 
are the ones that appear when the dimension of the space is taken into account. 
The essential reason for not including these relations is that applying properties 
of the building tensor do not generate this kind of relations. They arise when 
the whole monomial is taken into account so we need some tool to generate the 
dimensional properties for a given monomial. We are leaving this new tool for 
a future study. 
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Appendix: Mathematica Session. Demonstration 
of various properties relating Riemann monomi- 
als 

We present here an example of how the algorithm Simplif yAlllndex can be 
applied. We have chose the demonstration of well know polynomial properties 
of the Riemann tensor so the session will be also a short test for the algorithm. 

The session is part of a Mathematica notebook. The inputs are indicated by 
In [] : = . The outputs have no prefix. The titles and comments are enclosed 
using the symbols (* comment*). 

Since we use the Riemann tensor as example we introduce this tensor us- 
ing the TTC function InputSRiemann. This is internally equivalent to use 
InputTensor and InputSymmetries to declare the Riemann and Ricci tensors 
and to give the relation between them and between the scalar curvature. The 
enumeration of the various Riemann rules used in the session comes from [0 . 

On running the session we have used a Pentium-133 with 32Mb of RAM 
memory and the Mathematica version 2.2. The time to run the whole session is 
about 1500 seconds. The maximum memory used is about 2.6 Mb (the memory 
used to store the whole set of rules for future calculus is about 1Mb). The 
TTC program takes about 1.7Mb of memory so the total memory used by 
Mathematica is about 4.3Mb. 
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(*We have used the RiemannRule## list from the MathTensor book. 
In this way these is a good test*) 

In[] :=«ttc/ttc.m 



Tools of Tensor Calculus 4.0 
by A.Balfagon, P.Castellvi and X.Jaen 

\protect\vrule widthOpt\protect\href {http : / /balduf a . upc . es/ttc}{http : / /balduf a . upc . e 
e-mail: ttc@baldufa.upc.es 



TTCSettings= 
Compact [] =0f f 
Simplif yLevel=l 

TTCSimplify={SinTosin, Together, sinToSin} 

ScalarBasis=XX 

ScalarBasisQ=True 

Dimension [XX] =3 

$RecursionLimit=5000 

$SChristoff elTensor=G 



In[] := InputSRiemann [gn,XX, "R" ,Rie ,Ric ,R] 
{Rie, Ric, R} 

In[] := Inputlndex[{i, j ,k,l,m,n}] 
{i, j, k, 1, m, n} 
ORiemannRulel*) 

In[] := RRl=Rie[-a,a,-c,-d]//Index[gn] 
k 

+ R 
i j k i j 

In[] := RR1//Index[gn, Simplif yAlllndex] 


(*RiemannRule6*) 

In[] := RR6=Rie [-a,-b,-c,-d] Ric [a, b] //Index [gn] 
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k 1 

+ R R 
1 j k 1 i j 

In[] := RR6//Index[gn,SimplifyAllIndex] 


(*RiemannRule7*) 

In[]:= RR7=Rie [-a,-b,-c,-d] Rie[-e,a,b,c]+ 
1/2 Rie [-p,-d,-q,-r] Rie [-e ,p, q,r] //Index [gn] 
k 1 m 



R R 
j kilm klm 

+ + R R 

i j 2 j k 1 m i 



In[] := RR7//Index[gn,SimplifyAllIndex] 


ORiemannRulelO*) 

In[] := RR10=Ric [-a,-b] Rie [-c , -d, -e , a] Rie [b, c ,d, e] + 
1/2 Ric[-p,-q] Rie[-r,-s,-t,p]Rie[q,t,r,s]//Index[gn] 

j k 1 m i 



R R R 

jklm i ij lmk 

R R R + 

i j klm 2 



In[] := (RR10//Index[gn,SimplifyAllIndex]) //Timing 




(*RiemannRulel3 : here we have detected an error. In the Mathtensor 
book the second term has 1/4 as a factor. The correct factor is 
1/2*) 

In[] := RR13=Rie [-a,-b,-c,-d] Rie [-e , a, -g,b] Rie [c ,d, e ,g] - 
1/2 Rie [-p,-q,-r,-s] Rie [-t , -u,p , q] Rie [r, s ,t ,u] //Index [gn] 



k 1 m n i j 
R R R 
i j k 1 m n 

k 1 m n i j 

R R R 

i j k 1 m n 



12 



2 



In[] := RR13//Index[gn,SimplifyAllIndex] 




(*RiemannRulel5*) 

In[] := RR15=Rie [-a,-b,-c,-d] Rie[-e,a,-g,c] Rie [b,g,d,e] - 

Rie [-p,-q,-r ,-s] Rie [-t,p,-u,r] Rie [q,t , s ,u] + 

1/4 Rie [-p,-q,-r ,-s] Rie [-t , -u,p , q] Rie [r, s ,t ,u] //Index [gn] 

j m 1 n i k 
-(R R R ) + 

i j k 1 m n 

k 1 m n i j 

R R R 

i j k 1 m n 
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j m 1 n i k 
R R R 

i j k 1 n m 

In[] := RR15//Index[gn,SimplifyAllIndex] 


(*RiemaimRulel6*) 

In[]:= RR16=Ric[-a,-b, . ;b]-l/2 R[.;-a] //Index [gn] 
• ;i 

i R i . ; j 

+ R 

2 j 

In[] := RR16//Index[gn,SimplifyAllIndex] 


(*RiemannRulel8*) 

In[] := RR18=Ric[-a,-b, . ;-c, . ;b]- 



1/2 R[. ;-a, . ;-c] - Ric [-p, -a] Ric [-c ,p] +Ric [-p ,-q] Rie [-a,q,-c,p] //Index[gn] 
R 

. ;i . ; j k . ;k lk 

R R +R +R R 

ij 2 j ki ik.;j klij 
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In[] := RR18//Index[gn,SimplifyAllIndex] 


(*RiemannRulel9*) 
In[] := RR19=Ric [-a,-b, . ;-c, .; al- 
l/2 R[. ;-b, . ;-c]-Ric[-p,-b] Ric[-c,p]+ 
Ric [-p, -q] Rie [-b , q, -c ,p] //Index [gn] 



R 

. ;i . ; j k . ;k lk 

R R +R +R R 

ij 2 j ki ki.;j klij 



In[] := RR19//Index[gn,SimplifyAllIndex] 




(*RiemannRule24*) 
In[] := RR24=R[. ;-a, . ;-b, . ;a, . ;b] - 
1/2 R[. ;-p] R[. ;p]- 
R[. ;-p, . ;p, . ;-q, . ;q]- 
R[.;-p,.;-q] Ric [p,q] //Index [gn] 
• ;i 

-(R R ) 

• ;i • ;i • ; j • ;i • ; j 
R + R 

2 . ;i . ; j . ;i . ; j 

i j 

R R 

• ;i • ;j 

In[] := RR24//Index[gn,SimplifyAllIndex] 
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